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Universality in nontrivial continuum limits: a model calculation 
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We study numerically the continuum limit corresponding to the non-trivial fixed point of Dyson's 
hierarchical model. We discuss the possibility of using the critical amplitudes as input parameters. 
We determine numerically the leading and subleading critical amplitudes of the zero-momentum 
connected 2Z-point functions in the symmetric phase up to the 20-point function for randomly chosen 
local measures. Using these amplitudes, we construct quantities which are expected to be universal 
in the limit where very small log-periodic corrections are neglected: the [Z'-^''* (proportional to the 
connected 2Z-point functions) and the r2i (proportional to one-particle irreducible(lPI)). We show 
that these quantities are independent of the the local measure with at least 5 significant digits. 
We provide clear evidence for the asymptotic behavior U^-^'^* oc {2l)\ and reasonable evidence for 
r2t oc {2l)\. These results signal a finite radius of convergence for the generating functions. We 
provide numerical evidence for a linear growth for universal ratios of subleading amplitudes. We 
compare our r2t with existing estimates for other models. 

PACS numbers: 05.10.Cc,11.10.-z, 11. 10. Hi, 64.60. Ak 



I. INTRODUCTION 

Predictive theories are highly regarded because they 
offer many ways to be falsified. On the other hand, the- 
ories where every new observation results into the de- 
termination of a new parameter are uninteresting unless 
this process ends up with predictions. The condition of 
perturbative renormalizability allows us to get rid of UV 
regulators without a prohferation of adjustable parame- 
ters. It played a central role in the establishment of the 
standard model of electroweak and strong interactions. 
This model predicts many cross sections and decay rates 
Pi in terms of a very restricted set of input parameters 
and is considered as one of the most important accom- 
plishment of the 20-th century. 

In this article, we address the question of predictive- 
ness for the infinite cutoff limit associated with a non- 
trivial fixed point. This concept was developed with a 3 
dimensional example by Wilson In this case, we are 
free to use a much larger set of bare theories, including 
non-renormalizable interactions but we need to fine-tune 
one of the parameters in order to reach the critical hyper- 
surface, or stable manifold, which separates the ordered 
and disordered phases. The renormalization group (RG) 
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flows go along the stable manifold until they get close to 
the non-trivial fixed point, and then end up along the un- 
stable manifold. As the fined-tuned parameter gets close 
to its critical value, we can parametrize the zero mo- 
mentum (/-point functions in terms of power singularities 
multiplied by critical amplitudes. A few years ago 0, we 
have suggested to use (some of) the critical amplitudes 
as input parameters. However, it requires that we know 
how to pick an independent set of these amplitudes. 

Before the nineties 4] this question has not generated 
much interest. However, in the last decade, universal ra- 
tios of amplitudes associated with one-particle irreducible 
2Z-point functions at zero momentum have been calcu- 
lated with various methods 0, B Sll H IH [ill (see Ref. 
|12| for a review and a more extended list of references) . 
These ratios are denoted r2i as in Ref. where Table 
VIII summarizes the values of rg , rg and rio obtained in 
the literature. We are only aware of one calculation Q 
of ri2 and with large error bars. As we will see it 
is difficult to extrapolate the asymptotic behavior of the 
r2i from this data. 

In this article, we calculate universal quantities asso- 
ciated with the 2Z-point functions in a model where very 
accurate calculations of these quantities are possible up 
to the 20-point function. We use Dyson's hierarchical 
model IT^ which is very close to the 3 dimensional 
model used in Ref. 0. For definiteness, this model 
and our method of calculation are reviewed in section 
llTI In section IHII we discuss the numerical calculation of 
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the critical amplitudes for the zero-momentum 2Z-point 
functions. The definition of the renormalized couplings 
in terms of these amplitudes is discussed in section IIVI 
where we also explain why we expect some dimensionless 
quantities made out of these couplings to be approxi- 
mately universal |l5j | . In section we verify numer- 
ically that our expectations are correct and that in the 
infinite cutoff limit, dimensionless renormalized couplings 
^gj are in good approximation universal. The results 
presented here extend up to the 20-point function and 
allow us to study the asymptotic behavior. We found 
good evidence for a factorial growth comparable to what 
is found in Refs. [itIITsL [l9ll20l | for other models studied 
in the context of multiparticle production. 

In section IVII we show that the ratios of subleading 
amplitudes are in good numerical approximation uni- 
versal for the model considered here, as expected in 
general pit |2^ l23l| . In section IVIII we calculate the 
r2i, show that they are universal and of the same or- 
der of magnitude as those calculated in the literature 
for other models. Their asymptotic behavior is compat- 
ible with a factorial growth as suggested by "Grifhths 
analyticity" 5j . In the conclusions, we discuss the in- 
terpretation of these results and the applicability of the 
method to more realistic situations. 



II. THE MODEL 

Dyson's hierarchical model [TsLll^ l has a special kinetic 
term that keeps its original form after a RG transforma- 
tion. There is no wave-function renormalization in this 
model and rj—O. The RG transformation can be sum- 
marized by a simple integral formula that describes the 
change in the local measure after n RG steps W„(0), into 
the one after one more step: 

W„+i((/)) =C„+iexp((/3/2)02) 

X /d0V„((0-,^')ci/2)I^„((0-K0')ci/2) , (1) 



where C„+i is a normalization factor which can be fixed 
at our convenience. The model has one free parameter 
c — 2^~^/^ which can be adjusted in order to match 
the scaling of a free massless field in D dimensions. In 
the following we use the value c — 2^/^ corresponding to 
D = 3 exclusively. For this value of c, the non-trivial 
fixed point is known with great accuracy |24|. A list 
of accurate values of the eigenvalues of the linear RG 
transformation about this fixed point is given in Ref. [25j | . 
In the following, these eigenvalues are denoted Ai with 
Ai ~ 1.427 the largest and only relevant eigenvalue. We 
restrict our investigation to the symmetric phase of the 
model and use the accurate methods of calculations 
available in this phase. 

This integral formula is very similar to the approxi- 
mate recursion formula j27l | used in Ref. 01 . The main 
difference is that we integrate 2 field variables (keeping 



their sum constant) in one RG step instead of 2^ field 
variables. As a result, the change in the linear scale of 
the blocks is 2^/^ (instead of 2). This difference has no 
effect from a qualitative point of view. For a quantita- 
tive comparison of the two cases see Ref. '2^. A more 
explicit description of Dyson's hierarchical model and a 
more complete list of references can be found in Ref. [2^ . 

The RG transformation defines a flow in the space of 
local measures W{(j)) = exp(— V^(0)). We require parity 
invariance and V — > +oo when |0| —>■ +oo at a quadratic 
rate or faster. In the following calculations, the bare 
parameters will appear in a local measure of the Landau- 
Ginzburg (LG) form: 



(2) 



We have considered the (randomly chosen) possibilities 
given explicitly in Table IIVI provided in the Appendix. 
One of the choices is p ~ A and corresponds to non- 
renormalizable interactions in perturbation theory. In 
addition, we considered the Ising measure W{(j)) — (5((/)^ — 
1). For each local measure, we have varied the inverse 
temperature /3 in front of the kinetic term in order to 
reach a critical value Pc where a bifurcation is observed 
(see [3 for a more complete description). 

We have then studied numerically the flows of the lo- 
cal measures for values of /3 slightly smaller than their 
respective Pc- After n RG steps, the Fourier transform 
of the local measure 



Rn{k) = 1 + a„^i/c^ + an,2k'^ + 



(3) 



contains the information concerning the zero-momentum 
q-point functions. The connected parts are defined 
by 



with 



ln(i?„(fc)) =<,ifc^ + <,2fc + 



2^ sites 



We can then calculate 



{21) 



((^^2" sites 'P^) y 
2" 



(4) 



(5) 



(6) 



in terms of the a„^i. For instance, 

X(f) = 12(-a^,i + 2a„,2)(8/c2)", (7) 

and 

xi^^ = 240 (a^^i - 3 a„^i a„,2 + 3 a„,3)(32/c3)" . (8) 

For P < Pc the Xn have a well-defined limit when n 
becomes large that we call x''^'^ and that we now proceed 
to parametrize. 
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III. NUMERICAL DETERMINATION OF THE 
CRITICAL AMPLITUDE 

For values of /3 slightly smaller than (3c, we have 



+ A<^Jl cos (iu ln(/3e - /3) + 0^"'^) + 



,(9) 



with known parameters 72/, A (not to be confused with 
the gap exponent) and lu. We use the hyperscaling 
Hal values 



with hi 



and 



72/ = 7(5^ - 3)/2 



7 ~ 1.299140730159 



A ~ 0.42595 



(10) 

(11) 
(12) 



The possibility of log-periodic terms were first discussed 
in Refs. 0, They were identified in the high- 

temperature expansion |^ |32| of the model considered 
here. The frequency is 



27r 

h^ 



(13) 



(21) 

The amplitudes Aper. are however quite small, typically, 
they affect the 16-th significant digit of the susceptibility 
and it takes a special effort to resolve them numerically. 

As a starting point, we will calculate the sublead- 

(21) 

ing amplitude A\ . Using the parametrization (3 = 
fic — 10~^ and a small 5x, we obtain (neglecting the log- 
periodic corrections) 



A 



{21) [x'-'^^Kx + (5a;)10-(^+'^^)''2' - x^'^^\x)l{)~''^^'] 
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[10 



-A(5i 



(14) 



amplitudes. We follow the general procedure outlined by 
Wilson in Ref . ■ We consider a sequence L = 1,2... 
of models with (3 = (/3c — ^i^u) where u is positive but 
not too large. We introduce the increasing sequence of 
UV cutoffs 

Al = 2^Ar , (16) 

with Afl an arbitrary scale of reference. We follow the 
notations of Ref. 2] and L which is proportional to the 
logarithm of the UV cutoff should not be confused with 
the linear size of the system. We define the renormalized 
mass 



Ai 



x(2)(/3,-Ar^7.) 



Given that 



2^, 



(17) 



(18) 



the dependence on the UV cutoff disappear at leading 
order and one obtains 



(19) 



A 



(2) 



A^^\^{^)^ + LPC ^ 



with the log-periodic corrections 



LPC = Agl cos 



in u H — In — 
7 VAl. 

(20) 

In the infinite cut-off limit {L — > 00), the subleading cor- 
rections disappear. On the other hand, the LPC do not 
and we are in presence of a limi t cy cle with a cutoff de- 
pendence quite similar to Refs. [3J,|33- Strictly speak- 
ing the infinite cutoff limit does not exists, however, for 
practical purpose, the effects of the oscillations are so 
small that it introduces uncertainties that are smaller 
than the accuracy with which we establish the univer- 
sality. Consequently, these oscillations will be ignored in 
the following. 

We can now define U^'^\ the dimensionless coupling 
constants associated with the g-point function as 



Ar^u) 



qil+D/2)-D 



(21) 



As X becomes large enough, we observe a region where The constant of proportionality is fixed in such way that 



aIj^^ from 



A\ stabilizes. Once A ^ are obtained, it is easy to find if the conjecture [3g that (-l)'+i;^(20 > Q is correct, 

then t/t-^') > 0. We also introduce a power of /3 such that 
if we reabsorb f3 in the field definition, we get comparable 
quantities for measures with different (3c. In summary, for 
£1 = 3 



A 



(15) 



More details regarding the numerical methods can be 
found in Ref. 1331. 



IV. NON-PERTURBATIVE 
RENORMALIZATION 

In this section, we give a non-perturbative definition of 
the renormalized couplings. In the next section, we dis- 
cuss the relation between these couplings and the critical 



= ^_iy+l^(2l) (-^(2)^(3-5/)/2 ^3(1-0/2 ^ 1^22) 

In this definition, it is understood that x^^''' and x'^' are 
evaluated at the same B^_^ 

In section X of Ref. U5] , it is shown that when /? — > /3c 
(and the RG flows end up on the unstable manifold), the 
lj(2i) g^j.g universal in the approximation where the log- 
periodic oscillations are neglected. Conversely, if we find 
that in this limit, the J/^^'^ are approximately universal, 
it indicates that the log-periodic oscillations are small. 
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TABLE I: Universal values of f/'^'^* . 



21 uV^n* — 

4 1.505871 

6 18.10722 

8 579.970 

10 35653.8 

12 3.57769 lO*' 

14 5.31763 10* 

16 1.09720 10" 

18 3.00025 lO" 

20 1.04998 10^*^ 



V. UNIVERSAL COUPLINGS 

We now consider the infinite cutoff limit [Z^^'-** of [/(^') . 
In this hmit, the subleading corrections vanish. Due to 
the hyperscaling relation Eq. I|10(l . the leading singulari- 
ties cancel and we obtain 

j^(20* ^ (42))(3-50/2 ^3(1-0/2. (23) 

Using the previously calculated amplitudes we find that 
up to ^ = 10, the C/^^'^* are all positive and in good 
approximation universal. The values for particular mea- 
sures are given in the Appendix (Table IVII|) . The ap- 
proximately universal values are displayed in TableQlwith 
uncertainties of order 1 in the last digit. 

We have fitted InJ/'-'^-'* with a constant plus a linear 
term and a third term which is either \nq or hi[q\). Fig. 
^ shows two fits of these 9 values. The first fit (Fit 1 in 
Fig. HI) is 

;7(9'* ~ 21.5exp(3.436g)g""-^ , (24) 
and the second fit (Fit 2 in Fig. 

U^'J^* ~ 0.756(g!)i-^9 exp(-0.88g) (25) 

The two fits can barely be distinguished in Fig. If 
as in Fit 1 we use a general fit of the form 

~ A((j!)^ exp(Cg) , (26) 

then the values of the parameters change if we exclude 
the points with low values of q. For instance if we exclude 
the first five points, B ~ 1.15 instead of 1.29 if we use 
all the data. The intermediate values are shown on Fig. 
12 Using a nonlinear fit to determine how B depends on 
this choices, we found that B ~ 0.99 -I- O.Glg""'^^ where 
q is the smallest index in the dataset used to obtain B. 
We conclude that the asymptotic value is very close to 1 
and that the the leading growth is 

C/(«)* « q\ , (27) 

This result is similar to what is found in Ref . 0, ITsL [l9l 

I^J for other models studied in the context of multiparti- 
cle production. Note that the generating function of the 
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FIG. 1: ln([/''^'*) and the two fits described in the text. 
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FIG. 2: Value of the fitted exponent of the factorial for differ- 
ent sets of data. The abscissa is the smallest value of q used 
in the data set. 



connected 2Z-points function has a 1/{21)\ factor at order 
21 (see Eq. (jSJ) which indicates that the expansion of the 
generating function of the connected functions in powers 
of an external field has a finite radius of convergence (see 
section rVlII for more discussion). 



VI. CORRECTION TO SCALING AMPLITUDES 

Ratios of subleading amplitudes are also expected to 
be universal 0, 0| . To express these quantities we 
first define the relative strength of the corrections a^'' = 



5 



21 
4 
6 
8 
10 
12 
14 
16 
18 
20 



TABLE II: Universal values of S'^'** 




2.03 
3.26 
4.52 
5.8 
7.1 
8.3 
9.6 
11 
12 



wave function renormalization and consequently some of 
the rescalings introduced in CPRV are not necessary in 
the present case. 

In CPRV, a universal fmiction 



l>3 



(20! 



(28) 



obtained from the effective potential by suitable rescal- 
ings of the effective potential and its argument (the mag- 
netization) is defined (Eq. (5.7)). The universal coef- 
ficients r2i are ratios of amplitudes associated with the 
IPI. They can be expressed in terms of the connected 
Green's functions by expanding the magnetization as a 
power series of the magnetic field in the equation defining 
the Legendre transform. Explicit formulas for rg, rg and 
rio are given in Eqs. (5.23-25) of CPRV. These quanti- 
ties can be trivially reexpressed section For instance, 



re = 10 



([/(4)-)2 



(29) 



In CPRV, "Griffiths' s analyticity" is invoked to justify 
that A{z) in Eq. (|28|l has a finite radius of convergence. 
This requires, for large I, a growth not faster than 

r2i (X (Ci)'(20! . (30) 

The coefficients of r2i grow rapidly as can be observed 



FIG. 3: Ratios of subleading amplitudes S''^ 



om these, we define the ratios 

For the measures considered here, we found the approx- 
imately universal values shown in Table ^1 with uncer- 
tainties of order 1 in the last digit. Particular values are 
given in the Appendix. These values can be fitted well 
with a linear function as shown in Fig. 13 



VII. THE EFFECTIVE POTENTIAL 

In this section, we compare our results with related re- 
sults obtained for models in the same class of universality 
as the 3 dimensional Ising model with nearest neighbor 
interactions. In order to avoid repetitions, we follow sec- 
tions V, VII and Appendix C of Campostrino, Pelissetto, 
Rossi and Vicari 5j (CPRV for short). We follow exactly 
their notations except for the fact that we keep using 
superscripts in x'^'' to denote the connected 2g-points 
function (denoted X2q in CPRV). It should also be men- 
tioned that in the case of Dyson's model, there is no 



1401400 - 560560 



^(2)^(6) 



36036 



;^(2)2^(6)2 



17160- 



.(2)^^(8) 



X 



(4)- 



792 



(2)'^ v(6) v(8) 



X' ' x'- 'x 



220 



and 



X 

^{2)^^(12) 



(4)- 



(31) 



190590400- 95295200 



v(4)2 



10090080 



3203200- 



^(2)2^(6)2 
(4)4 



126126 



^(2)^^(6)^ 



X 

/(2)2^(8) 



X 



(4)t 



1716 



^(2)^^(8)^ 



360360 



X 



(2)-^ ^(6) ^(8) 



(4)'= 



50050 



X 

^^(2)^^(10) 



(4)^ 



2002- 



X 

,(2)4^(6)^(10) 



X 

;^(2)5^(14) 

Xi'if 



(4)t 



364 



^{2)^^(12) 



(32) 
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TABLE III: Universal values of r2i calculated numerically and 
predicted using the CPRV method for the value of p given in 
the last column. 



21 


r2i 


r2i predicted 


P 




6 


2.0149752 


1.81 


1.8946 




8 


2.679529 


2.47 


1.8946 




10 


-9.60118 


-10.1 


1.9218 


t 


12 


10.7681 


8.93 


1.9460 


u 


14 


763.062 


753 


1.9685 


a 


16 


-18380.8 


-1.76 X 10* 


2.3197 


H 


18 


1.5553 X 10^ 


1.62 X 10*^ 


2.1889 




20 


1.0374 X 10'^ 


1.04 X 10'^ 


2.1181 





The first (constant) and last (proportional to x'^'"*) terms 
of the r2i are expected to grow like (2/)!. The first term 



of the r2i (10, 280, 15400, 1401400, etc. 
can be obtained from the recursion 



denoted r, 



(0) 
21 ' 



m/3] 



.(0) 



2(i-m) (2;-3m)!m!6'' 



(33) 



B- 10 




4 8 12 16 20 

q 

FIG. 4: Values of In |r-g| versus q 

shown in Fig. ^ however the behavior is not as smooth. 
It is possible to obtain decent fits of the data of Fig. ^ 
with the parametric form 



with the initial conditions — r^^ — 1- A detailed 
analysis shows that this leads to a (2^)! growth (with 
power corrections). On the other hand, we also expect 
X*^^'^ ~ {2iy. in view of the numerical results of sect ion Ivl 
Obviously, if similar rates are found for the intermediate 
terms, then the bound of Eq. H30|l applies. 

The numerical values of ?'g,...r2o for the four mea- 
sures are given m the Table IVIIll in the appendix. The 
results show universality with the same kind of accuracy 
as in Section The universal values are summarized in 
Table Hill It should be noted that the numerical values 
of the r2i are orders of magnitudes smaller than some of 
the individual terms. For instance, for ri4, the sum of all 
the terms (753) is 5 orders of magnitudes smaller than 
the constant term (1.9 x lO'*). This requires minute can- 
cellations, and since from the discussion above, some of 
the terms grow as {21)1, this is a good indication that the 
bound of Eq. H3U|) should be satisfied. 

We have checked the consistency of our results by using 
the method proposed by CPRV to predict ?'2(;+ i), u sing 
r2,- ■ ■r2i as input. As one can see from Table PTll the 
two results are in reasonable agreement (relative errors 
of the order of 10 percent). One can also calculate r2(i+2) 
etc.. from the same input, but the agreement is not as 
good. It should noted that Eq. (7.21) in CPRV, for 
the intermediate parameter p, admits in general more 
than one positive root. The numbers given here have 
been obtained by using the smallest positive root also 
given in Table IIIII We also found that the equations 
leading to the prediction of rg and rg were identical. A 
detailed analysis shows that it is due to the fact that the 
magnetization exponent (usually denoted /3) is 7/4 |2^ 
for Dyson's model . 

The logarithm of \r2i\ is displayed in Fig. ^ The 
growth looks roughly similar to the growth of the L/'^^'-'* 



\n\rq\c^A + B\n{ql)+qC, (34) 

which is compatible with with Eq. (|30(l . However, 
the lack of smoothness makes the discrimination against 
other behavior difficult. 

A better way to identify the asymptotic behavior con- 
sists in studying the ratios 

p(2l) ^ r_2i+2 ^3j.^ 
r2i 

The factorial rate of Eq. (|30|l would imply that for I large 
enough, 

P(2') oc f , (36) 

which means a line of slope 2 in a log- log plot. Such a plot 
is provided in Fig. [S] The data is quite scattered, how- 
ever, the overall rate of growth seems compatible with 
Eq. I|36|) . the slope of the linear fit being 2.36. Note that 
without the last four data points, one might be tempted 
to conclude that the ratios are constant (this would imply 
an exponential growth rather than a factorial one). The 
same quantity with r2i replaced by C/(^')* is also provided 
for comparison. In this case, the linear behavior is quite 
evident and the slope slowly decreases from 2.22 to 2.11 
as we remove one by one the first five data points as was 
done in section IVl 

Our values of r2i are not very different from — 
2.048(5), rg = 2.28(8) and rio = -13(4) obtained by 
CPRV with an improved high-temperature method or 
the values obtained by other authors with other meth- 
ods (see Table VIII in CPRV for details). This can be 
explained from the fact that even though Dyson's model 
belongs to a different class of universality, the critical ex- 
ponents are not very different. This is encouraging for 
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a 




FIG. 5: Values of ln|p(20| 

versus ln2L The stars correspond 
to ratios of r2i as defined in Eq. (I35t . The hexagons corre- 
spond to the same quantity but with r2i replaced by JJ^^'^* . 
The lines are linear fits, the slope of the solid line is 2.11 while 
the slope of the dashed line is 2.36. 



the possibility of improving the hierarchical approxima- 
tion. Estimates of ri2 and rn can be obtained from Ref. 

Using the translation r2; — (2Z — l)!i^2i-i, we obtain 
ri2 = 38(16) with a sharp cutoff LPA and, ri2 = 20(12) 
and ri4 = 560(370) at lowest order in the derivative 
expansion. It is also interesting to compare the values 
Fs = 0.01679, F7 = 5.317 x lO"'' and Fg = -2.646 x 10"^ 
obtained with our data with the various entries of Table 
7 in Ref. Again, the values have the same order 

of magnitude, but no precise correspondence exists with 
any of the approximations listed. 



VIII. CONCLUSIONS 

We have calculated numerically the leading and sub- 
leading amplitudes corresponding to 4 randomly chosen 
local measures for the 2Z-point functions up to I —10. We 
found good evidence for approximate universal relations 
which allow, for the model considered here, to predict 
the amplitudes of all the connected 2/-point functions 
in terms of the amplitude for the 2-point function. We 
found clear indication that the universal amplitudes as- 
sociated with the connected 2/-point function grow as 
{21)1. In the context of perturbation theory, this growth 
seems to be related to the asymptotic nature of the ex- 
pansion fl7llT8LIT9ll20l |. However, our result is completely 
nonperturbative and indicates that the expansion of the 
generating function of the connected functions in power 
of the external field has a finite radius of convergence 
(since this generating function, as the universal function 
A{z) of Eq. (|2HIl, has 1/(2^)! factors in its definition). 

These results can also be used to calculate universal 
coefficients, denoted r2i, appearing in the effective po- 
tential and which have been calculated with a variety of 



methods H H S H H IH IH 113 for the universality 
class of the 3D Ising model with nearest neighbor in- 
teractions. Our results are compatible with universality 
with at least 5 significant digits. They are of the same 
order of magnitude as the existing estimates for other 
models. To the best of our knowledge, we have presented 
the first numerical estimates of rig, rig and r2Q. They 
are compatible with a factorial growth of the r2i and the 
expectation that the universal function A{z) has a finite 
radius of convergence j^. 

Our model does not pretend to be realistic. The fact 
that the non-perturbative continuum limit in this simpli- 
fied case leads to a situation where we have high pre- 
dictivity should be seen as an encouragement to look 
for nontrivial fixed point in more realistic theories. The 
application of the method to four dimensional models 
requires further investigation. The common wisdom is 
that in the pure scalar case, there is no non-trivial fixed 
point in 4D [38L Is^ (see however the discussion of Refs. 

El El)- On the other hand, for 4D models in- 
volving fermions and gauge fields, the question is more 
complex and stretches the limits of our present computa- 
tional abilities (see for instance the discussion of compact 
abelian fields coupled to fermions E^ or BCS inspired 
models of top condensation E3)- the Higgs sector 
will soon be probed at unexplored energies, a special ef- 
fort should be made to understand these questions. 
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APPENDIX: NUMERICAL RESULTS FOR 
PARTICULAR MEASURES 

In this appendix, we provide numerical results ob- 
tained with the measures given in Table Hvl and the Ising 
measure. 

TABLE IV: Measures used for particular calculations. 



LG(I) 


Exp[-O.5(;!)2-lO0*] 


/3c=7.7036412465997630 


LG(II) 


Exp[-2(?!)^-0.2(^'^] 


/?c=3.1273056619243551 


LG(ni) 


Exp[-0.1<^^-0.4<^*] 


/3c=2.2259466376795976 



TABLE V: a^'Va'^' Ratios for LG measures 





LG(I) 


LG(II) 


LG(III) 


q 














2 


0.10424 


1.00000 


0.21404 


1.00000 


0.31178 


1.00000 


4 


0.21156 


2.02961 


0.43401 


2.02776 


0.63339 


2.03152 


D 


0.33980 


3.25988 


0.69834 


3.26266 


1.10602 


3.25877 


8 


0.47075 


4.51615 


0.96793 


4.52219 


1.04700 


4.51281 


10 


0.60271 


5.78214 


1.24482 


5.81582 


1.80103 


5.77662 


12 


0.73514 


7.05257 


1.51927 


7.09805 


2.19654 


7.04517 


14 


0.86783 


8.32554 


1.79472 


8.38498 


2.59268 


8.31573 


16 


1.00068 


9.60004 


2.07087 


9.67519 


2.98897 


9.58678 


18 


1.13361 


10.8753 


2.34760 


10.9681 


3.38638 


10.8615 


20 


1.26664 


12.1515 


2.62528 


12.2654 


3.78335 


12.1347 



TABLE VI: a^"*' /a<^^ Ratios for the Ising measure 



<? 






q 






1 






2 


0.5614 


1.00000 


10 


3.24785 


5.78527 


18 


6.11352 


10.88978 


4 


1.13943 


2.02962 


12 


3.96229 


7.05787 


20 


6.83208 


12.16971 


6 


1.83040 


3.26042 


14 


4.66892 


8.31657 








8 


2.53632 


4.51785 


16 


5.39558 


9.61094 









TABLE VIL Dimensionless couplings (7^'''* 



q 


LG(I) 


LG(II) 


LG(III) 


Ising 


4 


1.5058710 


1.5058706 


1.5058710 


1.5058709 


6 


18.10722 


18.10721 


18.10722 


18.10722 


8 


579.9701 


579.9698 


579.9702 


579.9702 


10 


35653.80 


35653.77 


35653.80 


35653.80 


12 


3.577694E6 


3.577690E6 


3.577694E6 


3.577694E6 


14 


5.317628E8 


5.317622E8 


5.317628E8 


5.317627E8 


16 


1.097204E11 


1.097203E11 


1.097205E11 


1.097204E11 


18 


3.00025E13 


3.00024E13 


3.00024E13 


3.00025E13 


20 


1.04998E16 


1.04997E16 


1.04998E16 


1.04997E16 
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TABLE VIII: r2i values for the four measures considered, their averages and estimated errors. 



21 


LGI 


LGII 


LGIII 


Ising 


Average 


Error 


6 


2.0149752 


2.0149751 


2.0149752 


2.0149752 


2.01497516 


6x10"" 


8 


2.6795292 


2.6795279 


2.6795295 


2.6795289 


2.6795289 


7x10"'' 


10 


-9.6011836 


-9.6011888 


-9.6011824 


-9.6011849 


9.601185 


3x10"^ 


12 


10.768076 


10.768118 


10.768070 


10.768086 


10.76809 


0.00002 


14 


763.06239 


763.06185 


763.06229 


763.06192 


763.0621 


0.0003 


16 


-18380.758 


-18380.885 


-18380.933 


-18380.816 


-18380.85 


0.08 


18 


155520.69 


155541.95 


155500.45 


155543.79 


155526.72 


20.42 


20 


1.0366406E7 


1.0378039E7 


1.0377558E7 


1.0373698E7 


1.0374x10'' 


5x10^ 



10 



[1] K. Hagiwara et al. (Particle Data Group), Phys. Rev. 

D66, 010001 (2002). 
[2] K.Wilson, Phys. Rev. D 6, 419 (1972). 
[3] J. Godina, Y. Meurice, and M. Oktay, Phys. Rev. D 57, 

R6581 (1998). 

[4] V. Privman, P. C. Hohenberg, and A. Aharony, in Phase 
Transitions and Critical Phenonema, Vol. I4, edited by 
L. Domb and J. Lebowitz (Academic Press, New York, 
1991), pp. 4-121. 

[5] M. Campostrini, A. Polissetto, P. Rossi, and E. Vicari, 
Phys. Rev. E60, 3526 (1999), cond-niat/9905078. 

[6] N. Tetradis and C. Wetterich, Nucl. Phys. B422, 541 
(1994), hcp-ph/9308214. 

[7] R. Guida and J. Zinn- Justin, Nucl. Phys. B489, 626 

(1997) , hcp-th/9610223. 

[8] T. R. Morris, Nucl. Phys. B495, 477 (1997), hep- 
th/9612117. 

[9] A. Pelissetto and E. Vicari, Nucl. Phys. B522, 605 

(1998) , cond-mat/9801098. 

[10] M. Campostrini, A. Pelissetto, P. Rossi, and E. Vicari, 
Phys. Rev. B62, 5843 (2000), cond-mat/0001440. 

[11] M. M. Tsypiu, Phys. Rev. Lett. 73, 2015 (1994). 

[12] A. Pelissetto and E. Vicari, Phys. Rept. 368, 549 (2002), 
cond-mat/0012164. 

[13] F. Dyson, Comm. Math. Phys. 12, 91 (1969). 

[14] G. Baker, Phys. Rev. B 5, 2622 (1972). 

[15] Y. Meurice, Phys. Rev. E (in press) (2003), cond- 
mat/0312188. 

[16] G. Parisi, Statistical Field Theory (Addison Wesley, New 
York, 1988). 

[17] H. Goldberg, Phys. Lett. B246, 445 (1990). 
[18] J. M. Cornwall, Phys. Lett. B243, 271 (1990). 
[19] V. I. Zakharov, Phys. Rev. Lett. 67, 3650 (1991). 
[20] M. B. Voloshin, Nucl. Phys. B383, 233 (1992). 
[21] A. Aharony and G. Ahlers, Phys. Rev. Lett. 44, 782 
(1980). 

[22] M. Chang and A. Houghton, Phys. Rev. Lett. 44, 785 
(1980). 

[23] F. J. Wegner, in Phase Transitions and Critical Phe- 
nonema, Vol. 6, edited by L. Domb and M. S. Green 
(Academic Press, New York, 1976), pp. 7-124. 



[24] H. Koch and P. Wittwer, Math. Phys. Electr. Jour. 1, 

Paper 6 (1995). 
[25] J. Godina, Y. Meurice, and M. Oktay, Phys. Rev. D 59, 

096002 (1999). 

[26] J. Godina, Y. Meurice, M. Oktay, and S. Niermann, 
Phys. Rev. D 57, 6326 (1998). 

[27] K. Wilson, Phys. Rev. B. 4, 3185 (1971). 

[28] Y. Meurice and G. Ordaz, J. Phys. A (Letter to the Ed- 
itor) 29, L635 (1996). 

[29] J. J. Godina, Y. Meurice, and M. Oktay, Phys. Rev. D 
61, 114509 (2000). 

[30] T. Niemeijer and J. van Leeuwen, in Phase Transitions 
and Critical Phenomena, vol. 6, edited by C. Domb and 
M. Green (Academic Press, New York, 1976). 

[31] Y. Meurice, G. Ordaz, and V. G. J. Rodgers, Phys. Rev. 
Lett. 75, 4555 (1995). 

[32] Y. Meurice, S. Niermann, and G. Ordaz, J. Stat. Phys. 
87, 363 (1997). 

[33] M. B. Oktay, Nonperturbative methods for hierarchical 

models, (Ph. D. thesis), UMI-30-18602. 
[34] S. D. Glazck and K. G. Wilson, Phys. Rev. Lett. 89, 

230401 (2002), hep-th/0203088. 
[35] E. Braaten and H. W. Hammer, Phys. Rev. Lett. 91, 

102002 (2003), nucl-th/0303038. 
[36] J. Glimm and A. Jaffe, Quantum Physics (Springer- 

Verlag, New York, 1987). 
[37] T. R. Morris, Nucl. Phys. B458, 477 (1996), hep- 

th/9508017. 

[38] K. Wilson and J. Kogut, Phys. Rep. 12, 75 (1974). 
[39] M. Luscher and P. Weisz, Nucl. Phys. B290, 25 (1987). 
[40] K. Halpern and K. Huang, Phys. Rev. Lett. 74, 3526 

(1995) , hep-th/9406199. 

[41] K. Halpern and K. Huang, Phys. Rev. Lett. 77, 1659 

(1996) . 

[42] T. R. Morris, Phys. Rev. Lett. 77, 1658 (1996), hep- 
th/9601128. 

[43] J. B. Kogut and C. G. Strouthos, Phys. Rev. D67, 

034504 (2003), hep- lat/021 1024. 
[44] W. A. Bardeen, C. T. Hill, and M. Lindner, Phys. Rev. 

D41, 1647 (1990). 



